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A natural fueling mechanism that helps to maintain the main core deuterium and tritium (DT) 
density profiles in a tokamak fusion reactor is discussed. In H-moAe plasmas dominated by ion- 
temperature gradient (ITG) driven turbulence, cold DT ions near the edge will naturally pinch 
radially inward towards the core. This mechanism is due to the quasi-neutral heat flux dominated 
nature of ITG turbulence and still applies when trapped and passing kinetic electron effects are 
included. Fueling using shallow pellet injection or supersonic gas jets is augmented by an inward 
pinch of could DT fuel. The natural fueling mechanism is demonstrated using the three-dimensional 
toroidal electromagnetic gyrokinetic turbulence code GEM and is analyzed using quasilinear theory. 
Profiles similar to those used for conservative ITER transport modeling that have a completely flat 
density profile are examined and it is found that natural fueling actually reduces the linear growth 
rates and energy transport. 

PACS numbers: 28.52.Cx, 52.25.Fi, 52. 55. Fa, 52.65.-y 



The fueling of a tokamak fusion reactor is an important 
scientific and technological problem 0,1111. In ITER, 
deep central fueling is somewhat uncertain using present- 
day high-field side pellet injection schemes [J. Gas puff- 
ing and supersonic molecular beam injection [5j typically 
fuel the edge Q . In a Tokamak with deuterium, tritium 
and helium ash ion species, it was shown d that above 
certain level of density concentration, the helium ash nat- 
urally diffuses outward, helping in ash removal. Addi- 
tionally, due to quasi-neutrality, the main DT species go 
radially inward (or pinch) to balance the helium ash out- 
ward flux. In this paper, we show how adding cold fuel 
towards the edge can fuel the core and maintain the DT 
density profiles. Global multi-species gyrokinetic simula- 
tions using the GEM code 0|, show the effect that cold 
DT fuel added near the edge will naturally pinch towards 
the core. The main DT inward fluxes can be significantly 
reduced by fueling, while the outward flux of helium ash 
is maintained. The cold DT fuel will heat up due to ion- 
ion collisions as it moves towards the core, assuming the 
equilibration time is not too fast, the core density profiles 
may be sustainable with shallow pellet injection. We will 
show later that the transport time scale of the cold fuel 
is faster than the equilibration time in which it heats. 

Two different plasma profiles are studied here: the 
"Cyclone DIII-D base case" [1, ||| and a second case with 
ITER- relevant density and temperature profiles [II. A 
quasi-linear theory is presented to explain the direction 
of particle fluxes and the fueling mechanism. For the 
ITER-likc case, natural fueling has an additional benefit 
in that it is linearly stabilizing and reduces the energy 
transport of the main DT. We first present the simu- 
lation results for a Cyclone base case [1, [j|, which are 
parameters from a DIII-D if-mode assuming concentric 
circle flux-surface toroidal geometry. The "Cyclone base 
case" is a standard case widely used for benchmarking 
and for studying physics of turbulent transport in H- 
mode plasmas [9(. Here, we examine a global version 



of the Cyclone base case [8|], with kinetic electrons, and 
in the electrostatic limit. Recent studies [llj] show that 
increased collisionality can reduce the anomalous parti- 
cle pinch particle transport. We assume a collisionless 
plasma focusing on parameters towards the core region 
of a high temperature fusion reactor and the effects of 
collisionality will be investigated in a future publication. 

For the case presented, there are now five ion species: 
the main DT, the helium ash, and the DT fuel. The 
physical parameters are the major radius R = 445/9 s and 
minor radius a = 160p s , where p s = y/T e /TiPi and pi 
is the thermal proton gyroradius at a reference temper- 
ature T . The q profile is q(r) = 0.85 + 2.2 x (r/a) 2 . 
The density n(r) and temperature T(r) profiles, nor- 
malized to the reference density hq and temperature 
To, are shown in Fig. [T] The main DT profiles are se- 
lected so that L~}p(r) = B nt TSCch 2 [(r — i^/w], where 
L n = -n/(dn/dr), L T = -T/(dT/dr), r = a/2 and 
w = A0p s . The constants B n _T are determined by re- 
quiring R/\L n {r )\ = 2.2 and R/\L T (r )\ = 6.9. We 
assume the helium ash and electrons have the same 
temperature as the main DT. A linear density profile 
is assumed for the helium ash, as shown in the fig- 
ure. The density profile of the DT fuel has the form 
n(r) = C{l + tanh[5(r/a-0.65)]}/4, in which C denotes 
the density concentration, e.g., C = 0.05 means the peak 
density of the cold DT fuel is 5% of peak density of main 
core DT. This cold fuel is chosen so that the peak of 
the fuel density is towards the edge of the simulation do- 
main, at r/a = 0.8. It is important to note that this is a 
multi-component plasma and that the total DT profile is 
monotonically decreasing with minor radius throughout 
the simulation domain. 

The temperature of the DT fuel is set to be constant 
and is equal to the temperature of the main species at 
the simulation boundary near to the edge. Hence, the DT 
fuel is cold compared to the main DT species. The three 
dimensional (3D) simulation domain of 0.15 < r/a < 



2 



0.4 



0.3 



0.2 



0.1 



- main D/T 

- Helium ash 
fuel D/T 



0.2 



0.4 0.6 

r/a 




FIG. 1: (Color online). Density and temperature profiles of 
the Cyclone base case. Left: densities of main DT, helium 
ash and the DT fuel; right: temperature of the main species. 
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FIG. 2: (Color online). The Cyclone base case: the volume 
averaged particle fluxes of (a) the deuterium fuel, (b) the main 
deuterium, and (c) the helium ash for C = 10 -6 , 0.01, 0.02, 
0.03, 0.05, 0.07 and 0.1. The flux is normalized to nnVth, 
where Vth is the proton thermal velocity at To. 



0.85 is discretized using a 128x64 x 32 grid. The time step 
is At = 2f2~ 1 , where fii is the proton gyrofrequency at 
Tq. We use 4, 194, 304 particles per species with realistic 
mass ratios, e.g., the helium mass is 4 x 1837 times of 
the electron mass. The simulation uses GEM, which is a 
global gyrokinetic Sf particle-in-cell code @, EH ■ 

We study the effects of natural fueling by running sim- 



ulations with different concentrations of the DT fuel den- 
sity. Figure [2] shows the the particle fluxes of the deu- 
terium fuel, main deuterium, and the helium ash for con- 
centrations from C = 10 -6 (no fueling) to C = 0.1. The 
particle fluxes are calculated as T = ^ J2j WjVEx, where 
Wj is the weight of particle j, vex is the particle's E x B 
drift velocity in the r direction, and N is total number of 
particles. The fluxes are saved at 8 points along the mi- 
nor radius r, and the results shown in Fig. [5] arc volume 
averaged over all the points. The main and fueling tri- 
tium fluxes are similar to those of the deuterium for this 
case and not shown. From Fig. [2ja) , the cold fuel flow 
always goes inward, and the level increases for bigger C . 
Hence, "natural fueling," that is, the inward pinch of the 
cold fuel, is demonstrated for this simplified case. 
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FIG. 3: (Color online). The time-averaged particle fluxes of 
the main deuterium (MD), main tritium (MT), fuel deuterium 
(FD), fuel tritium (FT) and helium ash (HE) along r for C = 
0.01 and C = 0.05. 
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FIG. 4: (Color online) . The density-normalized particle fluxes 
Ta/rioVth of the five species for different C's. 

Estrada-Mila et al. @ predicts using both theory and 
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gyrokinctic simulation that for L n n e /L ne < 0.84, the he- 
lium ash flow should go outward, and this condition is 
well satisfied for the density profiles presented here. For 
the results in this paper, the helium ash always goes out. 
Assuming nearly adiabatic electrons, consistent with the 
properties of ITG turbulence, quasi-neutrality requires 
that the out-going helium ash flow must generate some 
inward going DT flow. Without the DT fuel, i.e., for the 
C = 1CP 6 case, the particle fluxes of main DT go inward 
as the helium ash goes out. This phenomenon is clearly 
shown in Fig. HJb) and (c). Therefore, edge fueling is 
needed to maintain the main DT density profile. 

By adding the cold fuel from the edge, the inward 
particle flux of the main DT is reduced. As shown in 
Fig.[2](b), for C = 0.05, the particle flux of the main deu- 
terium is nearly zero. Now, the helium ash goes out at 
the expense of the cold DT fuel instead of the hot main 
DT. Natural fueling is demonstrated and the main DT 
profile is maintained. The cold DT fuel, of course, must 
be heated as it migrates towards the core. 

This process is better demonstrated in Fig. [3j where 
the time-averaged particle fluxes of all five species are 
shown across the minor radius. Comparing the results 
of C = 0.05 to C = 0.01, the increased DT fuel fluxes 
apparently cancel that of the main DT. We note that al- 
though this cancellation is significant, it is not complete. 
This is because the positions of the maximum flux are 
different for the main and fuel DT, due to their individ- 
ual density profiles. The cold fuel concentration density 
should be kept below a threshold value above which the 
main DT flux becomes radially outward. As in Fig.^b), 
for the case C = 0.07 and C = 0.1, the main DT fluxes 
go outward in this case. 

Figure 0] shows the effects of fueling on all species. Af- 
ter normalizing by their individual densities, it is clear 
to see that fueling is most effective at pinching DT fuel 
inward. The particle flow of the main DT is altered from 
inward to near zero. The out-going helium ash flow ap- 
pears insensitive to fueling. This is important, in that 
the natural fueling does not increase helium ash build 
up. 

A rather crude assumption made in these simulations 
is that the cold fuel density is given and that it's temper- 
ature profile is constant. For this natural fueling mecha- 
nism to work, the fuel should remain cold as it migrates 
towards the core. The fuel temperature will come in equi- 
librium with the core temperature on the ion-ion equili- 
bration time. For a Tokamak fusion reactor at 15 keV in 
the core such as ITER, assuming the cold fuel temper- 
ature is 2 keV, the ion-ion thermal equilibration time is 
about 0.03 second. From Fig. [4] the averaged flow ve- 
locity of the fuel is about 7 x 10~ 4 «th, where Vth is the 
proton thermal velocity calculated at To. At fusion tem- 
perature, i>th is more than 10 6 m/s, therefore during the 
equilibration time, the fuel can go as far as 20m, which 
is much larger than the minor radius. The strong in- 



ward flux velocity is not simply caused by the negative 
density gradient of the fuel but rather due to the cold 
temperature. In fact, a pinch exists for cold fuel even in 
simulations with no density gradient. Indicating a con- 
vective pinch rather than diffusive density transport. 

Next, a multi-species quasi-linear theory Q is used to 
explain the natural inward pinch of the DT fuel. For 
simplicity, assume there are only two ion species: the 
hot main deuterium labeled as l i\ and the cold fuel deu- 
terium, as 'J'. The linear dispersion relation is Q 



Ri + Ri = 1 - iS, 



(1) 



in which the iS denotes the non-adiabatic electron re- 
sponse, 
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where kg = nq/r, n is the toroidal mode number, a = i, I. 
fi = e, fi = 1 — e, and e = nj/n e . The parameter 
u d = 2k e p s /R is a constant, u ma = k e p s /L nai u* t <t = 
keps/Lxa, and u* pa = cJ* n(T + ui*Ta- The quasi-linear 
particle flux given by 



r CT = Rc[ikg p s \cj)\ 2 R a {u)n e 



(3) 



is proportional to \<f)\ 2 at the nonlinear saturation level 
and the direction of the fuel flux is determined by 
— Im(i?/), with the inflow of the fuel corresponding to 
a negative Tj. 

The following conditions are generally satisfied in a 
tokamak plasma, and are assumed in this study: (i) 
Lu*Ti > 0, meaning the main ion temperature is higher 
towards the core, (ii) The theoryallows that for some 
Tokamak profiles, like in ITER [10j, u;*„, could be neg- 



ative. However, we require 



u t . 



U* T i > 0, 



so the temperature gradient should dominate, (iii) The 
fuel ion density should concentrate at the edge, therefore 
u*ni < 0. (iv) Since the fuel ions are cold, its temper- 
ature gradient should be small compared to its density 
gradient, and we require that u tp i = ui* n i + u>*ti < 0. 
Conditions (iii) and (iv) are essential for fueling. We can 
simply assume u)*ti = 0, as in simulations the fuel ion 
temperature is set to be the main ion temperature at the 
outer boundary point. 

Neglect the 0(kgp s ) 2 term and define u^ = (1 — 
+ eu mI , up = (1 — e)uj* pl + euj* p i, from Eq. (JTJ) we 
have Im(l/oj 2 ) = [(u N — w d )Im(l/w) + 5]/(uj d LUp). With 
Im(l/w) = — 7/|w| 2 where 7 is the linear growth rate, 
the — Im(i?/) term in Eq. (J3j) becomes 
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Since e is small, it is reasonable that up > 0. With 
the negative w* p /, the first term on the right hand side of 
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Eq. which is proportional to S, is therefore negative, 
meaning a non-adiabatic electron out flux should enhance 
fueling. The second term is proportional to uj^ which 
comes from the toroidal geometry of the device, and this 
term is also negative, so toroidal effects are also favorable 
for natural fueling. In the limit of lu*ti = 0, the third 
term becomes e(l — €)(jj m iuj*Til I (\w\ 2 u)p). Because of 
the density profile and the temperature gradient of the 
main ions, this term is also negative. We conclude that 
the fuel flux must be naturally negative and going inward 
towards the core. 

Due to quasi-ncutrality, adding a negative fueling flux 
should increase the main ion particle transport. This ef- 
fect can also be seen by switching I — > i and e — > 1 — e 
in Eq. j4j, in which all three terms on the right-hand- 
side are positive for the main ion, satisfying Im(i?j) + 
Im(i?/) = —5. In the case where only two ion species 
are present, e.g., hot main deuterium and cold fuel deu- 
terium, the fuel ion goes in at the expense of hot ions go- 
ing out, and therefore the natural fueling is of no value. 
However, when helium ash is present, the main DT pinch 
as the helium ash goes outward. Fueling is then useful to 
maintain the main DT profile, and balance the outgoing 
helium ash flux. 
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FIG. 5: (Color online). Profiles of the ITER-like case. 

The global version of the Cyclone base case presented 
here have monotonically decreasing density profiles that 
peak at the magnetic axis. So far there's no general con- 
sensus on the ITER density profile. One rather conserva- 
tive study [h]] shows ITER iJ-mode profiles with a com- 
pletely flat density. The inward peaking helium ash then 
causes the main DT density to be hollow, or peak near 
the edge. We now examine the natural fueling mecha- 
nism in this situation. The next set of simulations have 
the same physical parameters as in the Cyclone base case, 
but use the ITER-like hollow DT density and tempera- 
ture profiles, as shown in Fig. [5] 

The simulation results for this flat density ITER-like 
case are shown in Fig. [6] This completely flat density 
profile is overall much less unstable and as a result the 
cold fuel flux is also lower. This reduces the travel length 
of the fuel to about 0.8m before they are heated by ion- 
ion thermal equilibration. The outward going helium ash 
flow still causes the main DT to go inward without fu- 




FIG. 6: (Color online). The ITER-like case: the particle 
fluxes for C = 1(T 6 and 0.05. 



eling, but there is a different phenomenon. The main 
deuterium flux now separates from the main tritium flux. 
Similar to the previous case, adding cold DT fuel reduces 
the inward particle pinch of the main DT. However, due 
to the separation of the main D and T, it is very easy to 
overshoot with too large a C, thereby, causing the main 
tritium to go outward. This overshoot is shown in the 
right panel of Fig. [3J Given this effect of DT separation 
in the ITER-like case, one may consider using different 
density profiles for the deuterium fuel and tritium, possi- 
bly less tritium and more deuterium, but it is beyond the 
scope of this Letter and will require future investigations. 

Another significant new result of the ITER-like case 
presented here is that the linear growth rate of the 
ITG-induced instability is reduced by natural fueling, as 
shown in Fig. [7J This bonus is due to the hollow density 
profile of the main DT. In general, at low real frequency 
the linear growth rate 7 from Eq. fl} is roughly y/uJduJp 
with fueling and ^jujs^* V i for with main ions only. Since 
< lup < oj*pi, fueling decreases 7 and hence is a stabi- 
lizing effect. This effect is insignificant for the Cyclone 
base case: from C = ICP 6 to C = 0.1, 7 is only reduced 
by 2%, at about 6.8 x 10~ 4 fl,. However, in the ITER-like 
case, since io m i is negative, for comparable temperature 
profiles, l/u)*pi is bigger than the Cyclone base case and 



hence up/uj* 



1 — e + £Lu*pi/u}*pi is smaller. There- 



fore 7 is reduced more by fueling. As a result, the heat 
fluxes of both main DT and helium ash are also reduced. 
For C — 0.05 in the ITER-like case, the heat fluxes are 
reduced by nearly 50% from C = 10~ 6 . 
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FIG. 7: Linear growth rates for the ITER-like case. 
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